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ABSTRACT 

Context. Dynamical friction is expected to play an important role in a variety of astrophysical contexts, yet we still lack a quantitative 
understanding of this basic mechanism and of its effects. In fact, numerical simulations have shown that in inhomogeneous systems 
the classical idealized description, given by Chandrasekhar (1943), may have severe limitations, possibly related to the global nature 
of the processes that are involved. 

Aims. In this paper we address two issues. Firstly, we study how dynamical friction depends on the density concentration and on the 
pressure anisotropy of the host galaxy. For the purpose, we consider models characterized by a "realistic" distribution function and 
compare the behavior of dynamical friction in these systems to that found in other simpler models (often used in the past because of 
their mathematical convenience). Secondly, we study the response of the galaxy to the infall, by dynamical friction, of heavy objects 
("satellites") taken in a variety of initial configurations. 

Methods. The investigation is carried out by using a numerical laboratory set up in previous papers of this series. The process of 
dynamical friction is studied in terms of its strength (i.e., its efficiency to drag a satellite toward the galaxy center) and in terms of its 
ability to circularize the orbit of the satellite under friction. The response of the galaxy is studied in terms of the induced modifications 
to the galaxy density distribution and shape and of the changes produced to its phase space properties. 

Results. (1) We find that, within the range of our models, the pressure anisotropy present in the host galaxy has little effect on 
dynamical friction. Instead, the shape of the galaxy density profile is very important. The classical idealized description, although 
with an effectively smaller Coulomb logarithm, appears to be applicable to galaxy models characterized by a broad core (such as a 
polytrope) but not to concentrated models. Correspondingly, in contrast to the behavior found in models with a broad core, the orbits 
of satellites captured in concentrated models are not circularized by dynamical friction. To a large extent, these results confirm trends 
already known in the literature; in this respect, we also confirm the value of some simple modifications to the classical formulae 
proposed in the literature. (2) The induced evolution in the host galaxy reflects the initial conditions adopted for the captured satellite. 
Satellites spiraling in on quasi-circular orbits tend to modify the pressure tensor of the host galaxy in the tangential direction, while 
satellites captured along quasi-radial orbits tend to induce pressure anisotropy in the radial direction. While satellites captured along 
quasi-circular orbits make a galaxy change its shape from spherical to oblate, satellites captured along quasi-radial orbits tend to 
induce a shape in the host galaxy of the prolate type. This result suggests that the shape of early-type galaxies may just result from 
the characteristics of occasional mergers rather than being directly related to the effectiveness of the radial-orbit instability during the 
process of formation via collisionless collapse, as often argued in the past. 

Key words, methods: n-body simulations - galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: kinematics and 
dynamics 



1. Introduction 

Elliptical galaxies are essentially collisionless stellar systems. 
If they were purely collisionless isolated stellar systems, they 
might in principle remain in a dynamical equilibrium configu- 
ration for all their life. In practice, some evolution is expected 
to take place. In this series of papers (see Bertin et al. 2003, 
hereafter Paper I, and Arena et al. 2006, hereafter Paper II) we 
have investigated a process of dynamical evolution that should 
be ubiquitous. This is the evolution induced in the galaxy by the 
interaction with satellites or a population of heavy objects (such 
as a system of globular clusters or clumps of dark matter), which 
are dragged in toward the galaxy center by dynamical friction. 

The classical formula of dynamical friction was derived un- 
der highly idealised conditions (Chandrasekhar 1943). In con- 
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trast, real systems are inhomogeneous and generally charac- 
terized by a non-Maxwellian distribution function; in addition, 
satellites or other heavy objects dragged in by dynamical fric- 
tion are captured on complex orbits, on which resonant effects 
are expected. Therefore, in the last two decades one line of re- 
search has tried to determine, by means of semi-analytical tools 
or dedicated numerical simulations, to what extent the classical 
formula of dynamical friction is applicable to more complex and 
realistic situations. 

Theoretically, dynamical friction results from resonant ef- 
fects and is often interpreted in terms of interaction with a 
"wake" induced in the otherwise collisionless stellar system by 
the passage of the heavy object under investigation (e.g., see 
Mulder 1983, Tremaine & Weinberg 1984, Palmer & Papaloizou 
1985, Weinberg 1986, 1989; see also Sect. 2 of Paper I). An in- 
teresting connection has been drawn with the general statement 
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of the fluctuation-dissipation theorem (Kandrup 1983, Nelson & 
Tremaine 1999). Additional references on the subject have been 
recorded in Paper I and Paper II. 

With the help of self-consistent numerical simulations (see 
Bontekoe & van Albada 1987, hereafter BvA87, and Bontekoe 
1988, hereafter B88), it became clear that the process of dynam- 
ical friction depends on the model of the galaxy where it takes 
place. In particular, it was found that for a satellite dragged to- 
ward the center in a polytropic galaxy model (BvA87, confirmed 
in Paper I), the classical formula is basically applicable, although 
with a value of the Coulomb logarithm smaller than unity, and 
the orbits of satellites are circularized by dynamical friction; it 
was also noted that in these broad-core models dynamical fric- 
tion seems to disappear when the satellite reaches the central 
regions. In turn, in galaxies described by a King model (B88, 
Hashimoto et al. 2003) the classical formula fails (because one 
finds an effective Coulomb logarithm that varies with radius) and 
orbits are not circularized. 

It has been suggested that, if the maximum impact param- 
eter b max appearing in the definition of the Coulomb logarithm 
were chosen properly, a description in terms of the classical for- 
mula of dynamical friction might remain applicable. One choice 
(Hashimoto et al. 2003; Spinnato et al. 2003) has been b max = r s , 
where r s is the distance of the satellite from the center of the 
galaxy. Another proposed choice (Just & Penarrubia 2005) has 
been b max ~ p/\Vp\, so that the maximum impact parameter 
should be identified with the local scale of the density gradient. 
In practice, it is likely that the process of dynamical friction in 
inhomogeneous systems reflects some global effects that cannot 
be captured by variations on the classical formula. Therefore, it 
is important to widen the set of models on which the process of 
dynamical friction is studied by means of numerical simulations 
in order to see whether some general trends can be identified for 
the behavior of such a complex phenomenon. 

In general, past theoretical investigations have referred to 
objects moving through stellar systems characterized by an 
isotropic (quasi-)Maxwellian distribution function. In some 
cases (e.g., see van den Bosch et al. 1999, Jiang & Binney 2000, 
Spinnato et al. 2003, and Just & Penarrubia 2005), simulations 
are carried out on models constructed numerically by imposing 
a desired density profile and an assumed Maxwellian velocity 
distribution, evolved for a few dynamical times to reach approx- 
imate equilibrium. One notable exception, also because of its in- 
terest in the evolution of the host system, is a study of the process 
of dynamical friction on galaxies inside clusters (Binney 1977), 
considered as non-spherical systems supported by an anisotropic 
velocity distribution. 

Recently, a few papers have addressed the important issue 
of the orbital decay of objects treated as "live" stellar systems 
(Penarrubia et al. 2004; Fujii et al. 2005), with the possibility 
that the captured object is eventually disrupted by tidal interac- 
tions inside the host galaxy. These investigations certainly go in 
the direction of a more realistic study of the relevant astrophys- 
ical issues; yet, the complications and the number of parameters 
introduced make it even harder to extract from them a general 
statement on the ways dynamical friction proceeds. 

Most of these investigations have focused on the fate of the 
heavy object under friction, while in the present series of papers 
(see Paper I and Paper II) we have discussed the mechanism of 
dynamical friction in the more general context of the evolution 
of the stellar system where the processes occur. Recently, some 
cosmological issues have made this point of view even more im- 
portant. In fact, it has been noted that evolution induced by dy- 
namical friction (in contrast to the expectations from a scenario 



of adiabatic evolution; see Paper II) tend to lead to systems with 
softer density profiles (El-Zant et al. 2001; Paper I; El-Zant et 
al. 2004; Ma & Boylan-Kolchin 2004; Nipoti et al. 2004), with a 
significant impact on the current debate about the observational 
counterparts to the universal halo density profiles found in cos- 
mological simulations (Navarro, Frenk & White 1996; Moore 
et al. 1998; Ghigna et al. 2000; Navarro et al. 2004). This line 
of work has motivated the formulation of cosmological simula- 
tions based on the use of the Fokker-Planck approach (Ma & 
Bertschinger 2004; see also Evans & Collett 1997, Weinberg 
2001), for which it would be crucial to know whether dynam- 
ical friction can be, at least approximately, described in terms of 
local diffusion coefficients. 

In this paper we study the dependence of the process of 
dynamical friction (and of the response of the host galaxy) on 
the density and pressure anisotropy profile of the stellar system 
where dynamical friction takes place. The study thus goes in 
the direction of investigating the slow evolution of more real- 
istic models of elliptical galaxies and is intended to broaden the 
set of systems in which the process of dynamical friction, which 
is likely to have a global character, is studied quantitatively. To 
do so, we consider a family of stellar dynamical models, called 
/ (v) models, that has been widely studied, especially in relation 
to its ability to describe the products of collisionless collapse 
(Stiavelli & Bertin 1987, Bertin & Trenti 2003, Trenti & Bertin 
2005, Trenti, Bertin & van Albada 2005). This is basically a 
one-parameter family of models, with significant variations of 
density and pressure anisotropy profiles (from the concentrated 
models that reproduce the observed R 1/4 law to models with a 
broader core that are close to being unstable against the radial- 
orbit instability). Along this family of models, density concen- 
tration and pressure anisotropy change in a way determined by 
the physically justified distribution function. Therefore, in or- 
der to decouple the role of density concentration from that of 
pressure anisotropy, in this paper we also consider two series 
of models for which pressure anisotropy is specified and var- 
ied independently, by following the so-called Osipkov-Merritt 
prescription (Osipkov 1979, Merritt 1985): a Plummer model, 
as a prototypical broad core model, and a Jaffe (1983) model, 
as a prototypical concentrated model. For this study, we con- 
sider single satellites starting from a variety of initial conditions, 
and, following the arguments of Paper I and Paper II, shells of 
satellites that allow us to study the relevant mechanisms within 
a quasi-spherical environment. 

In Sect. 2 we present the models and the code. In Sect. 3 
we define the adopted units and the relevant diagnostics and we 
present some test simulations. In Sect. 4 we describe the results 
of the simulations for the process of dynamical friction in / (v) 
models compared with those observed in Plummer and Jaffe 
models with varying anisotropy. In Sect. 5 we give the results 
relative to the induced evolution in the previously considered 
galaxy models. In Sect. 6 we provide a general discussion and 
draw the main conclusions. 

2. The models and the code 

2.1. The host galaxy 

The basic models used in the numerical simulations to represent 
the initial state of the host stellar system are constructed from 
three different distribution functions. 

The first is the / (v) distribution function: 

/<"> = Ae~ {aE+d(fl millyll] , (1) 
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for E < and zero otherwise. The related models, constructed 
by imposing the self-consistent Poisson equation, have been 
briefly considered in Paper II and have been extensively de- 
scribed in a separate article (Trenti & Bertin 2005). Here, we 
recall that E and J are the single-star energy and angular mo- 
mentum and that the four real constants A,a,d, and v are pos- 
itive. A model is thus specified by two dimensional scales and 
two dimensionless parameters, such as v and the dimensionless 
central potential T* = -a4>(0). In the following, we take v = 3/4 
(the properties of the models would have only a mild depen- 
dence on v) and thus refer to a one-parameter family of models. 
All the models are characterized by a quasi-isotropic central re- 
gion and by radial pressure anisotropy in their outer parts. For 
higher values of ¥ the models are more concentrated and more 
isotropic, with properties similar to those of some observed ellip- 
ticals (in particular, they are characterized by an R 1 ^ 4 projected 
density profile) and of the products of collisionless collapse ob- 
tained in a variety of numerical experiments; for lower values 
of *P the models develop a relatively broad core and become 
more anisotropic, so that for *P < 4 they are unstable with re- 
spect to the radial-orbit instability (see Trenti & Bertin 2005 and 
Trenti, Bertin & van Albada 2005). In the following study, we 
would like to investigate the process of dynamical friction not 
only under "realistic" density profiles (higher values of *F), but 
also close to conditions of marginal stability of the host system 
CP* 4). 

The simulations that we will describe in Sect. [4] cover the 
range of models with *P running from 4.0 to 6.6. These mod- 
els have different density and anisotropy profiles. As mentioned 
above, as ¥ increases, the models are more concentrated and less 
anisotropic. However, given the finite size of the satellite and the 
fact that the change of the density profile with mostly affects 
the innermost region (see Trenti & Bertin 2005), all these mod- 
els will appear to have similar density concentration (which is 
higher than that characteristic of the polytropic model studied in 
Paper I and Paper II or that of the Plummer models) compara- 
ble to the density concentration of the Jaffe models (see below 
for further description). Therefore, simulations on models with 
varying ¥ will mostly test the behavior of dynamical friction 
while the pressure anisotropy of the host galaxy is varied. 

To study cases in which the pressure anisotropy profile can 
be varied independently of the density profile, we have then con- 
sidered the following two distributions functions. 

One is the fp distribution function (Merritt 1985): 



f P = A(-aQ) 



7/2 
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for - 1 /a < Q < and zero otherwise. The quantities E and J are 
here combined in Q — E + J 2 /(2r 2 ). This distribution function 
generates a Plummer density profile: 



3M 



An (b 2 + r 2 ) 5 ' 2 ' 
and an Osipkov-Merritt pressure anisotropy profile: 



a(r) = 2- 



<vg> + <j> 
<v?> 



2r 2 



2 ' 



(3) 



(4) 



where r a is the anisotropy radius defined as the radius where 
a(r a ) = 1 . The constant a is related to the gravitational constant 
G and the mass of the galaxy M by a = b/(GM) and A and b are 
two dimensional scales. The half-mass radius tm is related to b 



by b — V2 2 / 3 - \tm- With these models we will investigate the 
role of pressure anisotropy at fixed density profile, for the case 
of a broad-core distribution. 

The other distribution function (see Binney & Tremaine 
1987, Eq. (4-143), p. 241) generates a Jaffe (1983) density pro- 
file: 



P(r) 



M 



An r 2 (b + r) 2 



(5) 



and the above-described Osipkov-Merritt anisotropy profile. For 
this model b = Tm- With these models we will investigate the 
role of pressure anisotropy at fixed density profile, for the case 
of a concentrated distribution. Given the singularity present in 
the density distribution at r — 0, in our finite-resolution simu- 
lations the creation of an initial configuration according to the 
last model requires that we let evolve the initial configuration 
briefly, so that it settles into a quasi-equilibrium state close to 
the singular analytic model. 

2.2. Single satellites and shells of satellites 

We will study the orbital decay of single satellites and, as we did 
in Paper I and Paper II, of satellites distributed within a quasi- 
spherical shell. 

Each satellite is described by a rigid Plummer sphere, char- 
acterized by mass and radial scale M s and R s (i.e., we take 
M = M s and b = R s in Eq. (f3]l), respectively; we consider scales 
Mf and Rj when we refer to one element of a shell made of Nf 
fragments (in general we consider Mf — MJN/). 

We take single satellites moving initially on circular, quasi- 
circular, eccentric, and quasi-radial orbits. We discriminate 
among different levels of eccentricity of the initial orbit by re- 
ferring to the quantity Voj/v c , where vo., is the magnitude of the 
initial velocity of the satellite and v c is the velocity of the satel- 
lite that would correspond to a circular orbit at the chosen initial 
radius in the given galaxy potential; in any case, the velocity 
vector of the satellite at the beginning of the simulation has no 
radial component. 

We then consider two different configurations for "shells of 
fragments", both extensively described in Paper I and Paper II. 
The first configuration is defined by the following shell density 
profile: 



Psheiiir) = poexp[-4(r - r sheU (0)) /R 



shell* 



(6) 



for \r—r s heii(0)\ < R s heii and vanishing otherwise. Here r s h e ii(0) is 
the initial position of the shell where the density is po and R s h e ii is 
the shell half-thickness. In the presence of this shell density dis- 
tribution, it is necessary to recalculate the initial galaxy model 
as described in Subsection 3.2.2. of Paper I. The shell density 
profile is then sampled with Nf simulation particles (the "frag- 
ments") with mass Mf = M S heiilNf, where M s i w u is the total 
mass of the shell, initially placed on circular orbits. 

The second configuration is constructed by extracting the 
fragments directly from the distribution function of the galaxy 
following the distribution of (eccentric) orbits characteristic of 
the assumed distribution function (see Subsection 3.2.3. of Paper 
I). 

2.3. The code 

The galaxy evolution is computed with the same collisionless 
(mean field) particle-mesh code described in Paper II; for de- 
tails, see Trenti et al. (2005) and Trenti (2005). For the pur- 
pose of the present study, the choice of the code and the general 
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set-up of the simulations are justified in Paper I and Paper II. 
The galaxy distribution function is sampled with N particles by 
means of a Monte Carlo procedure. The same method is used to 
derive the position of the Nf fragments of the shell, for runs that 
involve the study of shells of satellites. Each simulation particle 
of the galaxy is subject to the mean field produced by the en- 
tire galaxy and to the direct action of the satellites present. Each 
satellite interacts directly with all the other satellites (when more 
than one satellite is present) and all the simulation particles rep- 
resenting the galaxy. 



3. Units, diagnostics, and test simulations 

3.1. Units 

The adopted units are 10 kpc for length, 10 n M o for mass, and 
10 8 yr for time. Thus, velocities are measured in units of 97.8 
km/s and the value of the gravitational constant G is w 4.497. 
We refer to models for which the total mass is M = 2 x 1O H M 
and the half-mass radius is = 3 kpc. Correspondingly, the 
dynamical time t d = GM 5/2 /(2K) 3/2 (h ere K is the total kinetic 
energy of the galaxy) falls in the range 0.18-0.25 10 s yr. 



3.2. Diagnostics 

The strength of dynamical friction is studied in terms of the 
effective Coulomb logarithm In A, as done in Paper I and in 
BvA87, and directly in terms of the friction coefficient y, de- 
fined by the relation ddf - -yv s , where ddf is the deceleration 
suffered by a satellite as a result of dynamical friction and v s 
is its velocity in the frame of reference locally comoving with 
the galaxy particles that are responsible for the dynamical fric- 
tion process. We will compare the values of these two quanti- 
ties, measured in the simulations, with those expected from the 
Chandrasekhar theory. 

The detailed steps for the computation of In A from the simu- 
lations were given in Sect. 5. 1 of Paper I. In particular, we recall 
that, following BvA87, the quantity In A is estimated from the 
relation 47rG 2 lnA = -v s (dE sat /dt)/(M 2 pcF(v s )), where E sat 
and M s are the energy and the mass of the satellite, pc represents 
the density of the galaxy, and F(v s ) is the fraction of particles 
with velocity, with respect to the local direction of the motion of 
the satellite, smaller than v s . The fraction F is computed on the 
basis of the actual distribution function, which is not assumed to 
be Maxwellian. 

As to the friction coefficient y, we may start from the basic 
equation for the rate of change of the satellite energy 



d£. 



dt 



= -yM s v s -v s ~ -yM s v 2 . 



(7) 



Therefore, extracting the "experimental" value of y from the 
simulations requires similar steps to those required for the ex- 
traction of In A and is actually simpler, because no sampling of 
the galaxy velocity space is needed. The friction coefficient is a 
dimensional quantity and thus, from our simulations, (1/y) will 
come out expressed in units of 10 8 yr. 

For a single satellite dragged in by dynamical friction on a 
quasi-circular orbit, when we will compare the above-defined 
measured values of In A and y to those expected from the 
Chandrasekhar theory, In Ac;, and ych, for simplicity we will re- 
fer to the values of the Coulomb logarithm and of the friction 
coefficient computed from the classical theory on the basis of 



the properties of the unperturbed galaxy at each radius r; sim- 
ilarly, for the velocity of the satellite entering the relevant for- 
mulae, we will refer to the circular velocity, as a function of r, 
expected in the adopted initial galaxy model from the relation 
v s (r) = y/rd® c /dr. 

In the classical theory the value of the maximum impact pa- 
rameter bmax, which appears in In Ach, has been subject to dif- 
ferent interpretations; furthermore, such theory was developed 
for point-like masses. In this paper, for point-like satellites, we 
take In Ac/, = In [3M/(2M S )], which coincides with the formula 
given by Chandrasekhar if evaluated at r = ru, with b max = 3rM- 
In turn, for extended objects, we refer to the formula suggested 
by White (1976), given that our satellites are characterized by a 
Plummer density profile: 



lnA c/l = 



1 



In 1 + 



R 2 



b 2 +R 2 

max s 



(8) 



where we set b max - 3r^; for simplicity, we keep the nota- 
tion In Ac;,, even though the latter formula was not provided by 
Chandrasekhar. [Note that the Coulomb logarithm for a point- 
like satellite depends on its mass, while for an extended object 
it depends on its radial size.] A satellite will be considered to be 
point-like when R s is smaller than the radius where the two ex- 
pressions of the Coulomb logarithm given above have the same 
value. 

The discrepancy between the measured value and the ex- 
pected value of the Coulomb logarithm will be often expressed 
in terms of the quantity A(r) = In A/ In Ac/,. Note also that in 
the classical theory the Coulomb logarithm and the coefficient 
of dynamical friction are related to each other through Eq. (0. 

As in Paper I, in our plots we will often use the dimension- 
less Lagrangian radial coordinate M(r)/M (based on the adopted 
initial galaxy model) instead of the radial coordinate r, so as to 
compare directly our results with those of BvA87. This choice 
has the effect of expanding the linear radial scale in the inner 
region, especially for radii smaller than the half mass radius r^- 

A global measure of the strength of dynamical friction is 
given by the fall time t/ a u of a satellite, or of a shell of satel- 
lites, relative to a given initial radius r s (t — 0) = ro, i.e., by the 
time taken to reach the center of the galaxy. 

We will also study whether orbits tend to be circularized by 
dynamical friction. The amount of circularization is calculated 
by checking the evolution of the ratio R m i„/R max of the pericenter 
to the apocenter of the satellite along its orbit. 

Finally, the effects of dynamical friction on the host galaxy 
are studied by following the evolution of several quantities: the 
galaxy density profile p(r), the pressure anisotropy profile a(r), 
the mean velocity profile (v)(r) (as is known, during the fall of 
a single satellite the galaxy may acquire some net rotation), the 
central density softening (Ap/p), = o, the shape parameters e and 
77, and the global anisotropy parameter k = 2K r jKj, where K r 
and Kj are the total kinetic energy of the galaxy in the radial and 
tangential directions respectively. The radial profiles are defined 
by averaging the relevant quantities inside spherical shells, given 
the fact that the galaxy remains quasi-spherical in the course of 
evolution. The definition of the pressure anisotropy profile a(r) 
extends that of Eq. (HJl to the case when mean motions, possi- 
bly induced in the galaxy, are present. The slight departure from 
spherical symmetry is then quantified by the two shape parame- 
ters, defined as e = bja and 77 = c/a, where a > b > c are the 
lengths of the axes of the inertia tensor (see Trenti, Bertin & van 
Albada (2005) and Trenti (2005)), so that e = 1 identifies the 
oblate case and 77 = e the prolate case. 
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Table 1. Test runs. 



Run 


Galaxy Model 


r a /r M 


k 


Tl 


yw ip _ 4 Q 


0.9 


1.84 


T2 


f M , x V = 4.6 


1.1 


1.70 


T3 


fiy) t m = 5.0 


1.2 


1.62 


T4 


= 5.4 


1.3 


1.57 


T5 


/ (r) ,T = 6.6 


1.7 


1.40 


T6 


Polytrope 


DO 


1.00 


T7 


Plummer 


1000 


1.00 


T8 


Plummer 


5 


1.06 


T9 


Plummer 


1.4 


1.43 


T10 


Jaffe 


1000 


1.00 


Til 


Jaffe 


5 


1.06 


T12 


Jaffe 


1.4 


1.32 


T13 


Jaffe 


0.75 


1.66 



Table 2. Equilibrium of the host stellar system. 



Run 


Ar/r|, 


Ar/r\ 2 


Ae/e 


A77/77 


Ak/k 


Tl 


0.14 


0.01 


0.16 


0.17 


0.04 


T2 


0.08 


0.01 


0.06 


0.06 


0.02 


T3 


0.05 


0.01 


0.016 


0.02 


0.01 


T4 


0.05 


0.01 


0.01 


0.02 


0.01 


T5 


0.04 


0.008 


0.007 


0.01 


0.01 


T6 


0.04 


0.007 


0.005 


0.005 


0.01 


T7 


0.03 


0.007 


0.005 


0.006 


0.008 


T8 


0.04 


0.009 


0.006 


0.008 


0.007 


T9 


0.03 


0.007 


0.006 


0.008 


0.01 


T10 


0.03 


0.01 


0.006 


0.007 


0.02 


Til 


0.04 


0.04 


0.006 


0.007 


0.02 


T12 


0.04 


0.05 


0.007 


0.008 


0.02 


T13 


0.05 


0.06 


0.006 


0.009 


0.02 



3.3. Test simulations 

In the following Section, we will show the results from a large set 
of simulations covering a wide range of physical conditions. In 
all these simulations and in the test runs described in the present 
subsection, the galaxy is sampled by N - 250, 000 simulation 
particles. Several tests performed earlier, and discussed in Paper 
I (see Sect. 5.1) and in Paper II (see Sect. 4), have demonstrated 
the adequacy of this number of simulation particles for the goals 
that we have set. 

Table[T]lists the properties of several test runs that have been 
performed to check the overall performance of the code in rela- 
tion to the models considered in this paper. In these test runs, a 
single satellite of very small mass M s = 10~ 9 M is placed ini- 
tially on a circular orbit at ro = 3.3/"^, with circular velocity v c 
(in physical units, this is in the range 260-300 km/s). The Table 
specifies the type of galaxy model that is used and the value of 
its total anisotropy parameter. 

The simulations are carried out for 30 dynamical times. 
Table [2] summarizes some observed limits on the galaxy evo- 
lution. The second column lists the maximum variation Ar/r|i 
of the radii of the spheres containing from 1% to 5% of the total 
mass of the galaxy; the third column gives the maximum vari- 
ation Ar/r\2 of the radii of the spheres containing from 50% to 
99% of the galaxy mass; the fourth and fifth columns list the rel- 
ative variations of the shape parameters and the last column the 
variation in the total pressure anisotropy. The significant varia- 
tions noted for run Tl reflect the fact that the equilibrium 
model with *P = 4 is close to the margin of the radial orbit insta- 
bility (Polyachenko & Shukhman 1981; Trenti & Bertin 2005). 

Table[3]summarizes the results of the test runs from the point 
of view of the stability of the satellite orbit and of the global 
conservations. The second column (Ar s /r s ) and the third column 
(AE sat /E sat) give the relative variations of the radial position and 
of the energy of the satellite; the last two columns represent the 
conservation of total energy and angular momentum per dynam- 
ical time for the entire system (galaxy plus satellite). 

At our disposal we also have a test simulation to check the 
equilibrium of the galaxy in the presence of a shell of fragments, 
i.e. run D3 of Paper II (with Nj = 25000 fragments of radius 
Rf = 0.33r M and mass Mf = 8 ■ 10~ 6 M); the results are in 
agreement with those of the corresponding model in the absence 
of the shell described above. 



Table 3. Stability of the orbit of a test particle and global con- 
servations. 



Run 


Ar s /r s 


^Esat i Asa/ 




AJiot 1 'J \ { 


Tl 


0.05 


0.02 


1.1 10" 3 


1.8 10- J 


T2 


0.07 


0.01 


8.4 10~ 6 


9.7 10" 5 


T3 


0.05 


0.01 


5.5 10- 6 


4.7 10~ 5 


T4 


0.06 


0.01 


3.6 10- 6 


5.3 10- 6 


T5 


0.06 


0.01 


5.4 10~ 6 


3.4 10- 5 


T6 


0.05 


0.01 


7.0 10~ 6 


4.3 10" 5 


T7 


0.05 


0.01 


4.1 10~ 6 


6.3 10- 6 


T8 


0.04 


0.01 


4.1 10- 6 


1.2 10- 5 


T9 


0.05 


0.01 


3.3 10- 6 


1.8 10" 5 


T10 


0.07 


0.02 


5.2 10- 6 


2.0 10" 5 


Til 


0.08 


0.03 


5.2 10~ 6 


1.2 10~ 4 


T12 


0.08 


0.03 


5.2 10- 6 


2.1 10- 5 


T13 


0.09 


0.03 


5.3 10- 6 


2.1 10- 4 



4. Dynamical friction in models characterized by 
different concentrations and various amounts of 
pressure anisotropy 

Table |4] lists the main characteristics of simulations of the fall 
of a single satellite in / (v) models, which are the main focus 
of interest of this paper. Similarly, Tab. [5] summarizes the prop- 
erties of other simulations performed on different models (see 
Sect. 12. Il l, for which the amount of pressure anisotropy can 
be varied independently of the density profile; for these runs, 
M s /M = 0.1 and ro/rM - 3.3. As indicated in the Tables, we 
have considered different combinations of mass M s , radius R s , 
initial position r , and initial velocity vo s of the satellite (relative 
to the galaxy mass M, half-mass radius tm, and circular velocity 
v c at r ). 

Table [6] lists the simulations for the case of shells of satel- 
lites (see Sect. 12. 2\ in / (v) models. Following the notation of 
Paper I, runs labeled by B refer to the first type of shell (with 
fragments on circular orbits) while runs labeled by BT refer to 
the second type of shell (with fragments on non-circular orbits 
extracted from the distribution function of the galaxy). The mod- 
els labeled by T for runs of type B are slightly different from the 
corresponding models used for simulations with a single satel- 
lite, because the equilibrium configuration has to be re-computed 
in the presence of the shell, as mentioned in Sect. 12.21 In all these 
runs the half-thickness of the shell is R s heii = 0.33^- 
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Table 4. Runs with a single satellite in selected fr v> models. 



Run 




MJM 




ro/r M 


VOs/Vc 


Purpose 


Fl 


5.0 


0.07 


0.33 


3.3 


0.96 


Varying 


F2 


5.0 


0.10 


0.33 


3.3 


0.96 


M s 


F3 


5.0 


0.14 


0.33 


3.3 


0.96 


and 


F4 


5.0 


0.07 


0.03 


3.3 


0.96 


R s 


F5 


5.0 


0.10 


0.03 


3.3 


0.96 




F6 


5.0 


0.14 


0.03 


3.3 


0.96 




F7 


5.0 


0.10 


0.06 


3.3 


0.96 




F8 


5.0 


0.10 


0.17 


3.3 


0.96 




F9 


5.0 


0.10 


0.06 


3.3 


0.99 




F10 


5.0 


0.005 


0.33 


1.0 


0.99 




Fll 


5.0 


0.10 


0.33 


1.0 


0.96 


Varying 


F12 


5.0 


0.10 


0.33 


2.0 


0.96 


ro 


F13 


5.0 


0.10 


0.33 


6.0 


0.96 




F14 


4.0 


0.10 


0.33 


3.3 


0.96 


Varying *P 


F15 


4.6 


0.10 


0.33 


3.3 


0.96 


(quasi 


F16 


5.4 


0.10 


0.33 


3.3 


0.96 


circular 


F17 


6.6 


0.10 


0.33 


3.3 


0.96 


orbit) 


F18 


5.0 


0.10 


0.33 


3.3 


0.99 




F19 


5.0 


0.10 


0.33 


3.3 


0.40 




F20 


4.0 


0.10 


0.33 


3.3 


0.50 


Varying T 


F21 


4.6 


0.10 


0.33 


3.3 


0.50 


and R s 


F22 


5.0 


0.10 


0.33 


3.3 


0.50 


(moderately 


F23 


5.0 


0.10 


0.17 


3.3 


0.50 


eccentric 


F24 


5.4 


0.10 


0.33 


3.3 


0.50 


orbit) 


F25 


6.6 


0.10 


0.33 


3.3 


0.50 




F26 


5.0 


0.10 


0.33 


3.3 


0.70 




F27 


4.0 


0.10 


0.33 


3.3 


0.30 


Varying *F 


F28 


5.0 


0.10 


0.33 


3.3 


0.30 


(eccentric 


F29 


6.6 


0.10 


0.33 


3.3 


0.30 


orbit) 


F30 


4.0 


0.10 


0.33 


3.3 


0.10 


Varying T 


F31 


5.0 


0.10 


0.33 


3.3 


0.10 


(quasi-radial 


F32 


6.6 


0.10 


0.33 


3.3 


0.10 


orbit) 



4. 7 . Local properties of dynamical friction 

In this subsection we focus on the mechanism of dynamical fric- 
tion, by measuring its local strength in the simulations by means 
of the coefficient of dynamical friction y and of the Coulomb 
logarithm In A, as defined in Sect. 13. 21 Some of the results, espe- 
cially the measured behavior in relation to the density concentra- 
tion of the host galaxy and a curious behavior noted with respect 
to the direction of motion of the satellite, will actually demon- 
strate that dynamical friction cannot be reduced to a purely local 
process as envisaged in the classical theory. 



Table 5. Runs with a single satellite in galaxy models with 
different combinations of density concentration and pressure 
anisotropy. 





Run 


Model r a /r M 


k 


RJr M 


V0 S /v f 




POl 


Pnl vtrnnp 




1 


0.33 


0.98 




P02 


Pnl vtrnnp 




1 


0.33 


0.1 




P03 


Polytrope 


CO 


1 


0.33 


0.5 




P04 


Polytrope 


CO 


1 


0.33 


0.7 




P05 


Pnl vtrrvnp 




1 


0.06 


0.98 




PL1 


Plnmmpr 

1 1L11JJ111L1 


1000 


1 


0.33 


0.98 




PL2 


Plnmmpr 


1000 


1 


0.33 


0.5 




PL3 


Plnmmpr 


5.0 


1.06 


0.33 


0.98 




PL4 


Plnmmpr 


5.0 


1.06 


0.33 


0.5 




PL5 


Plnmmpr 


5.0 


1.06 


0.33 


0.7 




PL6 


Plnmmpr 


5.0 


1.06 


0.06 


0.98 




PL7 


Plnmmpr 


1.4 


1.43 


0.33 


0.98 




PL8 


Plummer 


1.4 


1.43 


0.33 


0.5 




JA1 


Jaffe 


1000 


1 


0.33 


0.98 




JA2 


Jaffe 


1000 


1 


0.33 


0.5 




JA3 


Jaffe 


5.0 


1.06 


0.33 


0.98 




JA4 


Jaffe 


5.0 


1.06 


0.33 


0.5 




JA5 


Jaffe 


1.4 


1.32 


0.33 


0.98 




JA6 


Jaffe 


1.4 


1.32 


0.33 


0.5 




JA7 


Jaffe 


0.75 


1.66 


0.33 


0.98 




JA8 


Jaffe 


0.75 


1.66 


0.33 


0.5 




JA9 


Jaffe 


0.75 


1.66 


0.33 


0.7 




JA10 


Jaffe 


0.75 


1.66 


0.06 


0.98 


Table 6. Runs with a shell of N f fragments (R s heii 


jr M = 0.333). 


Run 


T 


N f M she „/M 


Rf/r M 


r s i,eii(0)/r M Purpose 


Bl 


5.0 


20 0.10 




0.33 


1.0 


Varying 


B2 


5.0 


100 0.10 




0.33 


1.0 


N f 


B3 


5.0 


20 0.10 




0.066 


1.0 


Varying 


B4 


5.0 


20 0.10 




0.033 


1.0 


R f 


B5 


5.0 


20 0.10 




0.017 


1.0 




B6 


5.0 


20 0.07 




0.033 


1.0 


Varying 


B7 


5.0 


20 0.14 




0.333 


1.0 


Rf and 


B8 


5.0 


20 0.14 




0.033 


1.0 




B9 


4.0 


20 0.10 




0.033 


1.0 


Varying 


B10 


6.6 


20 0.10 




0.033 


1.0 




Bll 


5.0 


20 0.10 




0.033 


0.5 


Varying 


B12 


5.0 


20 0.10 




0.033 


2.0 


r s hellW 


BT1 


5.0 


20 0.10 




0.33 


1.0 


Varying 


BT2 


4.0 


20 0.10 




0.33 


1.0 





4.1 .1 . The coefficient of dynamical friction 

In Fig. Q] the coefficient of dynamical friction measured in the 
simulations (left column) is compared to the expectations from 
the classical theory (right column), in a given / (v) model for the 
host galaxy with ¥ = 5. The top panels refer to an extended 
satellite (runs Fl - F3) of varying mass. The middle panels de- 
scribe the case of a point-like satellite (runs FA — F6). The bot- 
tom panels illustrate the case of a satellite with given mass, with 
varying radial size (runs F2, F 5, Fl, and F%). The observed pro- 
files of y(r) are systematically less steep than the corresponding 
profiles ych(r) expected from the classical theory and, in gen- 
eral, the observed dynamical friction is smaller than expected. 



On the other hand, the scaling of the coefficient of dynamical 
friction with satellite mass and radial size is in general agree- 
ment with that predicted by the classical theory. This confirms 
results already obtained in BvA87. 

In addition, no significant variations in the coefficient y(r) 
are observed by changing the model of the galaxy in the range 
= 4.0 - 6.6. This suggests that the coefficient of dynamical 
friction depends only very weakly on the pressure anisotropy 
content of the galaxy model, at least within realistic situations 
(as described by the anisotropy profiles a(r) characteristic of the 
family of / (y) models). 
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7 sim (r) in simulations 
R, 0.33r M M/M 
(extended) 




K,-0.033r M M/M 
(point— like) 




7 ( . h (r) expected by 

Chandrasekhar 
theory 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
M(r)/M M(r)/M 



Fig. 1. The coefficient of dynamical friction in the p v > galaxy 
model with *P = 5, for varying satellite mass and radial size. 
Observed coefficient of dynamical friction y(r) (left panels) 
measured in the simulations of the fall of a single satellite vs. 
expected coefficient jch(f) (right panels), as a function of the 
dimensionless Lagrangian radius M(r)/M. The observed coeffi- 
cient in the simulations is generally smaller than expected, with 
differences that become more significant in the inner regions. 
With respect to the classical theory, a similar scaling with satel- 
lite mass and radial size is noted. 



We have also checked that a single satellite and a fragment of 
equal mass and radius extracted from the spherical shell of run 
Bl feel the same amount of friction, at least for the radial range 
where such a comparison has been made. We recall that the study 
of dynamical friction for a shell of fragments has the advantage 
of offering a better controlled symmetry, but may include addi- 
tional collective effects related to satellite-satellite interactions. 

By considering the process of dynamical friction on eccen- 
tric orbits (runs F19 — F32) we have met a curious effect, il- 
lustrated in Fig. [2] Apparently, the dynamical friction felt when 
a satellite is "inbound", i.e. moving towards the center of the 
galaxy, is larger than the friction suffered when it is "outbound". 
In addition, the friction felt by the satellite when it falls towards 
the center increases with the eccentricity of the orbit. Therefore, 
dynamical friction appears to depend on the direction of the ve- 
locity vector of the satellite relative to the direction of the density 
gradient of the host galaxy. 

4.1.2. The Coulomb logarithm 

To compare the behavior observed in our simulations with other 
results reported in the literature, starting with BvA87, we have 
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Fig. 2. Coefficient of dynamical friction measured on eccentric 
orbits. The left panel refers to measurements recorded at instants 
when the satellite is inbound (i.e., moving towards the center of 
the galaxy) while the right panel records measurements made 
when the satellite is outbound (moving away from the galaxy 
center). The dotted line refers to the fall of a satellite along a 
quasi-circular orbit (run F18), crosses and dots to eccentric or- 
bits with initial velocity Vq s /v c of 0.7 (run F26) and 0.5 (runs 
F 22), respectively. 



0.5 



Si m illation 

Chand rasekhar 

Mcmii from she 1 ]; : ; I: 




0.2 



0.4 

M(r)/M 



0.6 



0.8 



Fig. 3. Coulomb logarithm measured in the p v > galaxy model 
with *F = 5. The dotted line represents A(r) measured in run 
F18, for a satellite initially placed on a circular orbit (a very 
similar behavior has been found in the corresponding runs for 
different / (v) models with varying V P). The dashed horizontal line 
indicates the expectations from the classical theory. The solid 
line shows the average value of A(r) found for the 20 fragments 
of simulation Bl. 



measured the strength of dynamical friction along the orbit of the 
satellite also in terms of the Coulomb logarithm. For / (v) galaxy 
models, the discrepancy between measured (In A) and predicted 
value (In Ac/,) depends on the location where the measurement 
is made, in the sense that A = In A/ In Ach = Mr), as illustrated 
in Fig. [3] 

Figure |4] as a counterpart to Fig. [2] illustrates the curious 
dependence, noted earlier in this section, of the friction strength 
on the direction of motion of the satellite for eccentric orbits. 
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Fig. 4. Coulomb logarithm measured on eccentric orbits. 
Symbols and overall layout are as in Fig. [2] Note that outbound 
satellites appear to feel a friction strength consistent with a con- 
stant In A, while inbound satellites feel a friction strength con- 
sistent with that characteristic of quasi-circular orbits. 



4.1 .3. Role of the density concentration and of the pressure 
anisotropy of the galaxy 

We have performed a similar analysis of the coefficient of dy- 
namical friction and of the Coulomb logarithm in other models, 
with the following results. 

The shape of the density distribution plays an important role 
in determining the strength of dynamical friction. In models with 
a broad core, such as the polytropic and Plummer models of runs 
PO and PL, dynamical friction is weaker than expected but can 
be reconciled with the classical theory by rescaling the friction 
coefficient or the Coulomb logarithm by a constant factor (see 
left panels of Fig. [5j. In turn, the study of the concentrated mod- 
els (runs J A) confirms the discrepancies already noted in the / (v) 
models (see right panels of Fig. O. 

Different amounts of pressure anisotropy at fixed density 
profile have negligible effects on the strength of dynamical fric- 
tion (see Fig.|5J. 

4.1 .4. Comparison with two modifications of the classical 
Coulomb logarithm proposed in the literature 

In Fig. [6] we compare the value of the Coulomb logarithm 
measured in our simulations with the prediction of two sim- 
ple modifications proposed by Hashimoto et al. (2003) and by 
Just & Penarrubia (2005). Adapted to the case of our satel- 
lites characterized by a Plummer density profile, the two for- 
mulae are, respectively, In A# = In [r/(1.4R s )] and In A jp = 
m [Gop('")/(l-3/?. s |Vp|)], where Qq ~ 1. For the two galaxy mod- 
els considered in this comparison, good agreement with the sec- 
ond formula is found by setting go = 2.72, while the first for- 
mula appears to fail in the inner regions of the galaxy. 

4.2. Global properties of dynamical friction 

In this subsection we focus on the mechanism of dynamical fric- 
tion, by measuring its global properties in terms of the fall time 
tfaii of a satellite, or of a shell of satellites, relative to a given 
initial radius r s (t — 0) = ro, and of the process of circularization 
of the satellite orbits. 



4.2.1. The fall time 

The fall time can be seen as the integral form of the coefficient 
of dynamical friction. Therefore, a consequence of the smaller 



Jaffe 




0.2 0.4 
M(r)/M 



Fig. 5. Coefficient of dynamical friction (top panels) and 
Coulomb logarithm (bottom panels) for models with different 
degrees of density concentration and pressure anisotropy. Left 
panels refer to the broad core Plummer models, with different 
contents of pressure anisotropy: isotropic (run PLl, solid line) 
and radially anisotropic with r a - lArM (run PLl, dashed line). 
Right panels refer to the concentrated Jaffe models: isotropic 
(run JAl, solid line) and radially anisotropic with r a = 0.15 tm, 
k = 1 .66 (run JAl, dotted line). The effect of pressure anisotropy 
on the strength of dynamical friction is negligible, in contrast 
with that of density concentration. 



values of y observed in the simulations (see Sect. 14.1.1) , with 
respect to the values predicted by the classical theory, is a greater 
fall time of the satellite. This is illustrated in the top panels of 

Fig.m 

The top left panel of Fig.|7]represents the orbital decay of a 
single satellite initially placed on a circular orbit as observed in 
our simulations (run F18; solid line) compared to the expecta- 
tions of the classical theory (dotted line); the observed fall time 
is about twice the expected value. The function rc/,(f) for the 
classical theory is computed using Eq. (0 with the classical ex- 
pression for ych, under the assumption that locally the satellite 
moves on a circular orbit. Here the host galaxy is described by 
an / (v) model with W = 5. 

In the top right panel of Fig. [7] the solid line represents the 
observed orbital decay of one of the 20 fragments of a spheri- 
cal shell for run Bl, in the same galaxy model. The dashed line 
represents the case of a single satellite of equal mass (M s = 
0.005M) and radial size (R s = 0.33r M ), for run F10. The dot- 
ted line gives the function rch(t) from the classical theory. Since 
the mass of the fragment (Bl) or of the satellite (FIQ) is 20 times 
smaller than that of the satellite illustrated in the left panel (F 1 8), 
the evolution observed is slower; note that a direct quantitative 
comparison between the two panels cannot be made easily, be- 
cause the initial location of the satellite is different in the two 
cases. Surprisingly, the fall of the fragment in the shell is signifi- 
cantly slower than that of a single satellite of the same mass (run 
F10) falling alone inside the galaxy. 
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Fig. 6. Comparison between the Coulomb logarithm measured in 
our simulations and the prediction of two formulae suggested in 
the literature. The left panel refers to the / (v) galaxy model with 
*P = 5 and the right panel to the Plummer model. Solid lines and 
long dashed lines represent the prediction by Hashimoto et al. 
(2003) and by Just & Penarrubia (2005), respectively. Dotted 
lines represent lnA(r) measured in runs F18 (left) and PL1 
(right), for a satellite initially placed on a circular orbit. Short 
dashed horizontal lines indicate the expectations from the clas- 
sical theory. 



The bottom left panel of Fig.Qrepresents the orbital decay of 
a single satellite (run F22) initially placed on an eccentric orbit 
(vos = 0.5v c ). The bottom right panel represents the fall of a 
shell (run BTl), extracted from the galaxy distribution function, 
with fragments on a variety of eccentric orbits. The solid lines 
identify the radii of the spheres containing (from the bottom) 15, 
35, 55, 75, and 95 % of the mass of the fragments, respectively. 

Figure [8] illustrates the dependence of the fall time on var- 
ious properties of the satellite and of the host galaxy. The two 
top panels show the dependence on satellite mass for extended 
(left) and point-like (right) satellites. The middle left panel il- 
lustrates the dependence on the radial size of the satellite. The 
bottom panel shows the dependence on the eccentricity of the 
satellite orbit. As shown in the middle right panel, no significant 
differences are found by varying the value of T* of the / (v) model 
considered. In the top four panels, filled triangles refer to the ex- 
pectations of the classical theory and dots to the values observed 
in the simulations. 

In a series of polytropes with increasing size of the central 
core, B88 had noted that a single satellite is actually unable to 
reach the center in the course of its orbit. We found a similar 
behavior in simulations of models with a broad core (both in 
a polytropic model, P05, and in a Plummer model, PL6). For 
single satellites in f^ models this effect is not observed, prob- 
ably because the size of the central core in these concentrated 
models is too small with respect to that of the satellite; on the 
other hand, in simulations of shells of fragments initially placed 
on circular orbits in concentrated / (v) models, the shell does not 
actually reach the center of the galaxy, but rather settles down 
into a quasi-equilibrium configuration of finite size (see bottom 
right panel of Fig. |7|. Possibly, this latter behavior results from 
the fact that in the final stages the shell of fragments collectively 
mimicks the presence of a broad core. 

4.2.2. Circularization of orbits? 

Dynamical friction may circularize initially eccentric orbits (see 
Tremaine, Ostriker & Spitzer 1975). This effect has been ob- 
served in some simulations of broad core galaxies (e.g., in a 
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Fig. 7. Orbital decay produced by dynamical friction. Top left: 
fall of a single satellite (run F18, solid line) compared to the ex- 
pectation from the classical theory (dotted line). Top right: fall 
of a fragment in a shell (run Bl, solid line) and of a single satel- 
lite of same mass and radius (run F10, dashed line), compared 
to the expectation from the classical theory (dotted line). Bottom 
left: fall of a single satellite initially placed on an eccentric orbit 
(run F22). Bottom right: fall of a spherical shell of fragments 
(run BTl; the solid lines from the bottom represent the radii of 
the spheres containing 15, 35, 55, 75, and 95 % of the mass of 
the shell). 



polytrope, BvA87), but has been shown to be absent in the more 
concentrated King models (see B88). To contribute to the study 
of this problem, we have performed several numerical experi- 
ments in which the satellite starts from the outer regions of a 
galaxy described by a concentrated / (v) model, with a velocity 
vector in the tangential direction, but with speed smaller than 
that of the corresponding circular orbit at the same initial posi- 
tion. We have then followed the evolution of the ratio R m i„/R ma .x, 
computed along the orbit covered by the satellite. This study 
shows that, in practice, within the family of / (v) models orbits 
are not circularized by dynamical friction. 

One example is given in the top panel of Fig. |9j where the 
evolution of the ratio R m i„/R mnx is shown for the orbit of the 
single satellite of run F22, starting with vq s /v c = 0.5 and falling 
in an / (v) model with *P = 5. The orbit initially increases and then 
slightly decreases its eccentricity, in practice with no evidence 
for circularization. 

Changing the properties of the host galaxy within the fam- 
ily of / (v) models (runs F20, FIX, F22, F24 and F25), at fixed 
initial eccentricity, does not lead to significant changes in the 
observed circularization. In turn, we note that, by starting from 
quasi-radial orbits, some circularization appears to take place 
(compare the result for run F31 with those for runs F2, F22, 
F26, and F28 in the middle panel of Fig. 

We have also looked for the process of circularization dur- 
ing the fall of a spherical shell made of Nf = 20 fragments on 
eccentric orbits (run BTl). The results are illustrated in the bot- 
tom panel of Fig. [9] which displays the distribution of the values 
of Rmm/Rmax for the 20 fragments in the simulation. The quan- 
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Fig. 8. Dependence of the fall time on various properties of the 
satellite and of the galaxy. All panels display the fall time of a 
single satellite starting at j*o = 3r«- The top two panels represent 
the variation of t/ a n on the satellite mass for an extended (left 
panel, runs Fl, F2, and F3) and a point-like (right panel, runs 
FA, F5, and F6) satellite. The middle left, the middle right, and 
the bottom left panels show, respectively, the dependence of the 
fall time on the radial size of the satellite (runs F 2, F5, Fl, and 
F8), on the parameter *P characterizing the / (v) galaxy model 
(runs F2, F14, F15, F16, and F17), and on the eccentricity of 
the sateUite orbit (runs F2, F18, F22, F26, F2S, and F31). In 
the top four panels, filled circles refer to results from the simu- 
lations, triangles to the expectations from classical theory. 



tity iV(> Rmin/Rmax) on the y-axis is the number of fragments 
with ratio R m i„/R max greater than the value shown on the x-axis. 
The simulation was stopped at time f/,„ = 120^. By comparing 
the initial distribution of points (filled squares) with that of the 
final configuration (crosses), we can conclude that the number 
of fragments with eccentric orbits actually increases during the 
simulation. 



4.2.3. Role of the density concentration and of the pressure 
anisotropy of the galaxy 

The study of additional models with different density profiles 
and different amounts of pressure anisotropy basically confirms 
the general trends noted so far. 

For single satellites initially placed on quasi-circular or- 
bits, the fall time in radially anisotropic models is somewhat 
shorter than in models with the same density distribution but 
with isotropic pressure. In particular, the fall time in the most 
anisotropic Plummer model (r a = 1 ArM', more anisotropic mod- 
els are unstable) is 92% of that in the isotropic model; in the 
most anisotropic Jaffe model (r a = 0.75^) it is 78% of that 
in the corresponding isotropic model. The fall time in less con- 
centrated models is slightly longer than in more concentrated 
models. 

In broad core models (polytrope and Plummer), initially ec- 
centric orbits are almost fully circularized by dynamical friction, 
independently of the amount of pressure anisotropy present; in- 
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Fig. 9. Circularization of orbits in the model with *F = 5. Top 
panel: evolution of the orbit "aspect ratio" R m i„/R ma . x during the 
fall of a single satellite (run F22). The x-axis counts the number 
of turns made by the satellite starting from t — 0. Middle panel: 
initial (filled squares) and final (crosses) data points for the case 
of the fall of a single satellite on eccentric orbits (runs F2, F22, 
F26, F28, and F31). Bottom panel: initial (filled squares) and 
final (crosses) distribution of the orbit aspect ratio for a shell of 
satellites on eccentric orbits extracted from the distribution func- 
tion of the galaxy (run BT1); here A^(> Rmin/Rmax) is the number 
of fragments with orbit aspect ratio greater than the value given 
on the jc-axis. 



stead, in concentrated models (Jaffe) circularization does not 
take place (see Fig.fTOt. as observed in / (v) models. 

We may expect that the orbit of a satellite will be circular- 
ized if the rate of loss of angular momentum is negligible with 
respect to the rate of loss of energy; in the opposite situation, the 
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Fig. 10. Circularization of orbits in different galaxy models. 
Independently of pressure anisotropy, orbits are circularized in 
broad core models (runs PL2 and PL8) and remain eccentric in 
concentrated models (runs JA2 and JA6). The x-axis counts the 
number of turns made by the satellite starting from t = 0. 



orbit is expected to become more eccentric. As to the relative 
loss rates, in the simulations we have observed a rather complex 
behavior, as represented in the E s -J s plane. Here E s = E sa ,/M s 
is the satellite specific energy and J s is its specific angular mo- 
mentum. Relating diagrams of this type to some simple theoret- 
ical expectations (based on the structure of orbits for the unper- 
turbed spherical potential associated with the galaxy) is difficult, 
because in the simulations the satellite has finite mass and the 
galaxy is evolving. Figure QT| thus compares the orbits of two 
runs for a single satellite falling in a broad core model (Plummer; 
left frame) or in a concentrated model (/ (v) with — 5; right 
frame). In each galaxy model, the two runs represent a case of 
an initial quasi-circular orbit (vo s /v c « 1; dotted line) and a case 
of an initially eccentric orbit (Voj/v c = 0.5; solid line). The or- 
bits in the E s -J s plane are followed by the satellite moving from 
the upper right to the lower left part of the plane. The effect of 
circularization in one case and of lack of circularization in the 
other case are thus demonstrated. 



5. The evolution of the host galaxy induced by 
dynamical friction 

As a result of the interactions that determine dynamical fric- 
tion on a satellite or on a shell of fragments, the galaxy slowly 
evolves. 

The models of the / (v) family that we have considered range 
from systems with relatively high *P, with projected density pro- 
files well fitted by the R 1 * 4 law characteristic of the luminosity 
profiles of elliptical galaxies, to systems with *P = 4. For smaller 
values of T the models would be unstable with respect to the 
radial orbit instability (Polyachenko & Shukhman 1981; Trenti 
& Bertin 2005). One reason why we have decided to study the 
problem of dynamical friction for models close to the margin of 
the radial orbit instability is that we would like to test whether, 
under such conditions, the response of the galaxy to the presence 
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Fig. 11. Orbits in the specific energy - specific angular mo- 
mentum plane. The decay of an initially quasi circular orbit 
(vo s /v c » 1 ; dotten line) is compared to that of an initially eccen- 
tric orbit (vo. 5 /vc = 0.5; solid line), for a broad core galaxy model 
(left frame) and for a concentrated model (right frame). Along 
each curve, the satellite moves from the upper right to the lower 
left of the plane. The specific satellite energy E s = E sat /M s and 
angular momentum J s are given in code units. 



of infalling satellites may be significantly enhanced by internal 
collective effects. 

We have found that a single satellite and a shell of frag- 
ments induce significantly different forms of evolution in the 
host galaxy. In addition, we have found the unexpected result 
that a single satellite, under appropriate conditions, can signifi- 
cantly alter the state also of firmly stable systems. In all this, we 
should keep in mind that in the simulations of the fall of a single 
satellite, in general we are considering a rather heavy satellite 
(M s = 0.1M); clearly, lighter satellites are expected to produce 
less prominent effects. 



5.1. Evolution induced by the fall of a single satellite on firmly 
stable systems 

The effects of the fall of a single satellite on the density profile of 
the host galaxy are approximately independent of the eccentric- 
ity of its orbit (an example is given in Fig. [12] where the dashed 
line refers to a quasi-circular orbit and the solid line to a highly 
eccentric orbit). The evolution of the density profile of the host 
galaxy is in the direction of a softening of the initial density con- 
centration (see also discussion in Paper II). 

On the other hand, within the firmly stable part QV > 4) of 
the sequence of / (v) models, the final state attained by the galaxy 
in the case of the capture of a single satellite does depend signif- 
icantly on the initial eccentricity, as measured by vq s /v c , of the 
satellite orbit. In fact, in contrast with the quasi-independence of 
the final density profile from the eccentricity of the orbit of the 
satellite, the effects on the pressure anisotropy profile are sizable 
(see the example given in Fig. [13), and so are those on the final 
shape of the galaxy. 

Therefore, we have identified two different types of behav- 
ior that split the *F - vo s /v c plane in two regions separated by a 
transition boundary. In the first region (that we will call region of 
"negative feedback" evolution), the galaxy evolves slowly only 
due to dynamical friction, while in the second region (that we 
will call region of "positive feedback" evolution) evolution ap- 
pears to be governed by the combined effect of dynamical fric- 
tion and of the radial orbit instability. A schematic description of 
the available regimes is given in Fig. [14] 
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Fig. 12. Evolution of the density profile of the host galaxy in- 
duced by the fall of a single satellite. The initial galaxy density 
profile, multiplied by the volume factor so as to better illustrate 
the mass distribution, is given by the dotted line. The final den- 
sity profiles of the galaxy, after the fall of a single satellite on a 
circular (dashed line, run fT8) and on a highly eccentric (solid 
line, run F3 1 ) orbit, are similar. 



5.1.1. The "negative feedback" region 

The evolution of the galaxy in the negative feedback region is il- 
lustrated in the right column of Fig. [15] which uses the examples 
offered by the behavior of the / (v) model with T* = 5. 

The capture of a satellite on a quasi-circular orbit leads to 
final configurations characterized by a central isotropic core sur- 
rounded by a radially anisotropic envelope (see top right panel 
of the Figure). The global anisotropy parameter k decreases to- 
wards a final value smaller than the initial one. The final config- 
uration of the galaxy is characterized by rotation that in the inner 
parts (out to m r^) is rigid (as observed in the polytropic model; 
see Paper I), and in the outer parts is differential and tends to 
disappear (see bottom right panel). During evolution, the shape 
of the galaxy changes into an oblate configuration (middle right 
panel). 

Note that the net effect of the satellite is to "sweep away" 
the radial anisotropy present, and thus goes in the direction of 
removing a possible source of radial orbit instability. The final 
shape generated basically reflects the transfer of angular momen- 
tum from the satellite to the galaxy. 

5.1.2. The "positive feedback" region 

In this region the satellite is captured on highly eccentric or 
quasi-radial orbits (see left column of Fig. [TBI). 

The final pressure anisotropy profile differs substantially 
from that attained in the negative feedack regime. The pressure 
anisotropy profile remains basically unchanged, with a slight in- 
crease of radial anisotropy in the central region (top left panel). 
Indeed, the global anisotropy parameter k increases. The trend is 
thus in the direction of making the radial orbit instability active. 
The final shape is that of a prolate spheroid (middle left panel) 
and this can be seen as a combined effect of the transfer of energy 
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Fig. 13. Evolution of the pressure anisotropy profile induced on 
the p v > model with Y = 6.6 by the fall of a single satellite on 
orbits characterized by different eccentricities. The initial pres- 
sure anisotropy profile is given by the dotted line. The solid line 
represents the final profile for a capture on a highly eccentric or- 
bit (run F32); the dashed line shows the final profile for a fall 
on a quasi-circular orbit (run F17). In the latter case, note the 
development of an inner highly isotropic core surrounded by a 
tangentially biased shell, inside the external radially anisotropic 
envelope. 



and momentum from the satellite to the galaxy and the possible 
excitation of the radial orbit instability. Little or no rotation is 
noted in the final configuration (bottom left panel). 

5.1.3. The transition region 

In this region the heavy object moves on intermediate eccentric- 
ity orbits. The qualitative behavior of the final density profile is 
similar to that found for other conditions (see Fig.lT2T>. The pres- 
sure anisotropy profile remains basically unchanged (top central 
panel of Fig. Q3J the global anisotropy parameter remains con- 
stant). The system reaches a final configuration with a generic 
triaxial shape (middle central panel) and little or no rotation (bot- 
tom panel). 

5.1.4. Varying ¥ 

Higher *P models are more isotropic and more concentrated and 
tend to be changed less in the course of the evolution from their 
initial state induced by the interactions with a satellite. 

A curious behavior occurs in high-T* models. As we noted, a 
single falling satellite on a quasi-circular orbit, tends to change 
the pressure anisotropy in the tangential direction. On such 
high-*!* models, which have only small amounts of pressure 
anisotropy in the radial direction, the effect of a satellite is 
such that evolution can lead to a configuration characterized 
by the presence of a tangentially biased anisotropic shell posi- 
tioned between the inner isotropic core and the external radially 
anisotropic envelope (see Fig.[T3]l. 

In the negative feedback regime, the flattening of the oblate 
shape induced in the galaxy becomes smaller as one moves to 
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Fig. 15. Evolution induced by a single satellite in the various regimes identified in the "P — Vq s /v c plane. The left column (run F31) 
illustrates evolution in the region of positive feedback. In this regime, the pressure anisotropy profile is kept almost unchanged (top 
panel), while the galaxy evolves into a prolate spheroid (middle panel); some systematic motions in the radial and azimuthal direc- 
tions are detected in the final quasi-equilibrium configuration, but covered by noise (bottom panel). The central column describes 
the evolution in the transition regime (run F28), in which the galaxy reaches a generic triaxial shape. The right column (run F2) 
describes evolution in the region of negative feedback; in this regime, the the galaxy develops a quasi-isotropic central core, an 
oblate shape, and significant systematic rotation in the azimuthal direction. This figure illustrates a cut of Fig.[14]at *T = 5. 



models that are characterized by higher values of *P; correspond- 
ingly, in the positive feedback regime the prolate spheroid that is 
generated becomes less elongated (see Fig. [Toll. In other words, 
models away from the margin of the radial orbit instability are 
generally "harder", i.e. more resistant to changes, while models 
closer to *P = 4 are "softer", i.e. more vulnerable. 



5.1 .5. Role of the density concentration and of the pressure 
anisotropy of the galaxy 

Independently of the presence of pressure anisotropy, the evo- 
lution of the density profile induced by a satellite on a quasi- 
circular orbit is qualitatively similar for all models considered in 
this paper (polytrope, Plummer, Jaffe, and / (v) ); quantitatively, 
models with increasing density concentration are affected less. 
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Fig. 14. Different regimes of evolution induced by the fall of a 
single satellite in p y > galaxy models. The diagram refers to a 
satellite of radial size R s = 0.33?m and mass M s - 0.1M. In the 
upper region ("negative feedback") the galaxy evolves slowly 
towards an oblate configuration. In the lower region ("posi- 
tive feedback") the galaxy changes its shape significantly into 
a prolate spheroid. In the transition region the galaxy reaches a 
generic triaxial shape. 



For a given initial density profile, the final density profile of 
more anisotropic models is slightly shallower (runs PL1, PL3, 
PL1, JAl, J A3, JA5, JA1). 

When a satellite is dragged in on a quasi-circular orbit, all 
models reach a final oblate shape. For a given initial density pro- 
file, the oblate product is flatter if the model is more anisotropic 
to begin with. 

These results can be summarized by stating that more con- 
centrated and more isotropic stellar systems are affected less by 
evolution induced by dynamical friction. These results confirm 
the behavior showed by / (v) models. 

5.2. Evolution induced by a spherical shell of fragments on 
firmly stable systems 

The fall of a spherical shell (with fragments on circular or on 
eccentric orbits) induces a slow evolution in the density and in 
the pressure anisotropy profiles of the galaxy similar to that ob- 
served in the negative feedback region of the ¥ - vo. 5 /v e plane 
in the case of the fall of a single satellite. However, at variance 
with that case, the final configuration remains quasi-spherical, 
non-rotating, and characterized by a smoother anisotropy pro- 
file. One example of evolution of the anisotropy profile for the 
/ (v ' model with T* — 5 is given in Fig.fTTl 

5.3. Evolution for models at the margin of the radial orbit 
instability 

The / (v) model with *P = 4.0 is at the boundary between stable 
and unstable models (with respect to the radial orbit instability). 
Its evolution, as a result of interactions with a satellite or a shell 
of satellites, basically follows the same trends noted for firmly 
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Fig. 16. Final configurations induced in some f^ galaxy models 
by a satellite falling on orbits of different eccentricity. The top 
panel displays the ratio between the final and the initial global 
anisotropy parameter defined in Sect. 13.21 the middle panel the 
final value of e (e = 1 corresponds to an oblate shape), and the 
bottom panel the ratio between the final values of e and rj (e/r/ = 
1 corresponds to a prolate shape). Points on the left part of each 
panel are in the positive feedback region, those on the right are 
in the negative feedback region. Dashed and long-dashed lines 
refer to the final configuration induced by a spherical shell of 
eccentric satellites on the *T = 5.0 (run BTl) and ¥ = 4.0 (run 
BT2) models respectively. 
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Fig. 17. Evolution of the pressure anisotropy profile induced on 
the galaxy by the fall of a spherical shell of fragments. The solid 
line refers to evolution induced by a shell of fragments on quasi- 
circular orbits (run Bl) while the dashed line to evolution by a 
shell of fragments on eccentric orbits (run BTl). The final pro- 
files are very similar in the two cases. 
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Fig. 18. Evolution of the pressure anisotropy profile induced on 
the most radially anisotropic model simulated (*¥ — 4.0). Fall 
on a quasi-radial orbit (run F3Q) leads to some increase of the 
pressure anisotropy in the radial direction inside the half-mass 
radius and a corresponding decrease in the external region. Fall 
on a circular orbit (run F14), instead, leads to the formation of 
a quasi-isotropic core of radius just smaller than Tm, surrounded 
by a radially anisotropic envelope. 



stable models in the previous section, with no evidence of sud- 
den or discontinuous response. 

In particular, the fall of a single satellite on a quasi-circular 
or lowly eccentric orbit can "sweep away" the radial anisotropy 
present in the central regions and thus stabilize the galaxy 
against the radial orbit instability. Correspondingly, the total 
amount of anisotropy k decreases and the shape of the galaxy 
becomes oblate, with some systematic rotation. The pressure 
anisotropy profile becomes characterized by the particular shape 
visible in Fig. [18] where there is a sharp transition, at about the 
half-mass radius, between the central isotropic core and the ex- 
ternal radial envelope. This shape recalls that found in some 
equilibria studied by Trenti & Bertin (2006), which is associ- 
ated with models with global content of anisotropy k above the 
usually accepted threshold (Polyachenko & Shukhman 1981) for 
the onset of the radial orbit instability. 

Instead, satellites on highly eccentric and quasi-radial orbits 
bring the galaxy into conditions of instability (see filled circles in 
Fig.[l6]for vo.s/vc ~ 0.1). The value of k is almost unchanged (top 
panel), but the content of radial anisotropy in the inner region 
increases and the shape of the galaxy becomes prolate (bottom 
panel). 

In contrast, during the fall of a quasi-spherical shell of satel- 
lites, the galaxy maintains its round shape and absence of rota- 
tion, while evolving to a less concentrated and more isotropic 
configuration, as for more stable models (long-dashed line 
in Fig. [TBI). 

6. Discussion and conclusions 

In this paper, by means of N-body simulations, we have ad- 
dressed the problem of dynamical friction in a realistic envi- 
ronment (the / <v) models), i.e., in galaxy models characterized 



by significant density gradients and significant anisotropy in the 
velocity distribution. To understand the relative role of the two 
factors, density concentration and pressure anisotropy, indepen- 
dently of each other, we have also studied additional models with 
tunable (Osipkov-Merritt) pressure anisotropy profiles, for the 
cases of a Plummer and of a Jaffe density distribution. The prop- 
erties of dynamical friction have been analyzed in terms of the 
effects on the orbit of a satellite (or a shell of fragments) and of 
the corresponding evolution induced in the stellar system. 

We have found that the density concentration of the host 
galaxy has a significant impact on the strength of dynamical fric- 
tion, while the presence of pressure anisotropy appears to be less 
important. In broad core models, the process of dynamical fric- 
tion can be described by the classical theory of Chandrasekhar, at 
least approximately, with a smaller value of the Coulomb loga- 
rithm. In contrast, concentrated models, better suited to describe 
the density distribution of some real elliptical galaxies, are diffi- 
cult to reconcile with the classical theory, because the effective 
value of the Coulomb logarithm would formally change with ra- 
dius. Fortunately, this behavior appears to be reasonably well 
reproduced by the use of a simple modification of the classical 
Coulomb logarithm, proposed by Just & Penarrubia (2005). In 
this respect, and in view of current discussions of cusp forma- 
tion and evolution in self-gravitating systems, it should be noted 
that the variation of the empirical Coulomb logarithm in the in- 
nermost regions (M(r)/M < 0.2) turns out to be rather modest, 
even for concentrated systems. In addition, while in broad-core 
models, dynamical friction tends to circularize the orbit of cap- 
tured satellites, in concentrated models it does not. One curi- 
ous finding, irrespective of the density concentration of the host 
galaxy, is that apparently the dynamical friction felt by satellites 
falling on eccentric orbits is stronger for inbound than for out- 
bound satellites. 

As to the evolution induced in the galaxy by the fall of heavy 
objects by dynamical friction, we have found the following re- 
sults. Dynamical friction leads to a decrease in the central den- 
sity of the host galaxy. The fall of a single satellite makes a 
galaxy evolve into different final shapes and phase space proper- 
ties. In particular, a satellite on a quasi-radial orbit induces evo- 
lution towards a prolate shape and to more (radially) anisotropic 
configurations, while a satellite on a quasi-circular orbit leads 
to an oblate shape with pressure anisotropy changed in the tan- 
gential direction and with some rotation. The induced soften- 
ing in the density distribution and the changes in the pressure 
anisotropy are stronger in the case of a single satellite and for 
less concentrated and more radially anisotropic systems; they are 
less pronounced in more concentrated and more isotropic mod- 
els. 

The experiments performed in this paper show that the shape 
of an otherwise collisionless stellar system can be significantly 
modified by the capture of a single satellite of finite but relatively 
small mass, with the final configuration being oblate, prolate or 
even triaxial, depending on the orbital characteristics of the en- 
counter. Therefore, we argue that this process (which, in certain 
regimes, may couple to the stability properties of the galaxy with 
respect to the radial orbit instability) should be considered as one 
important cause for the distribution of shapes among elliptical 
galaxies, while, traditionally, such distribution was ascribed to 
the role of instabilities in more abstract terms, especially in terms 
of processes taking place during formation via collisionless col- 
lapse (e.g., see Aguilar & Merritt 1990, Cannizzo & Hollister 
1992, Theis & Spurzem 1999, Warren et al. 1992). We have 
checked that the softening of the central density profile, induced 
by dynamical friction, does not affect substantially the projected 
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density profile of the host stellar system that remains well de- 
scribed by the R 1/4 law; this agreement is better in more concen- 
trated models, because they are originally closer to the R 1/4 law. 
In other words, this study of dynamical friction may have led 
to identifying one other important role of minor mergers in de- 
termining the evolution of galaxies into their currently observed 
morphologies. 
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